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Abstract 

We investigate the way big rips are approached in a fully inhomogeneous description of the 
space-time geometry. If the pressure and energy densities are connected by a (supernegative) 
barotropic index, the spatial gradients and the anisotropic expansion decay as the big rip is 
approached. This behaviour is contrasted with the usual big-bang singularities. A similar 
analysis is performed in the case of sudden (quiescent) singularities and it is argued that the 
spatial gradients may well be non-negligible in the vicinity of pressure singularities. 
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Consider the situation where the matter content of the present Universe is dominated by 
a perfect fluid with barotropic index w — p/p < —1. Evidence of this possibility seems to 
be suggested from the analysis of Type la supernovae. If this is the case future singularities 
may be expected [1, 2] (see also Refs. [3, 4, 5]). 

In the present paper we intend to study the nature of future big rip singularities in a 
fully inhomogeneous approach whose relevance in the context of usual big-bang singularities 
has been exploited long ago [6] (see also [7, 8] and references therein). For instance, in the 
case of conventional big-bang singularities one can show that the relative contribution of the 
gradients decays as the singularity is approached. Does the same happens in the case of 
future (big-rip) singularities? In the case of big-bang singularity the anisotropy is believed 
to play an important role in the way the singularity is effectively approached, is this true 
also for big-rips? These are some of the questions we ought to address. 

Consider first the case where the perfect barotropic fluid filling the Universe is charac- 
terized by a supernegative equation of state, i.e. w = — 1 — e, with e > 0. We shall then 
be interested in the contribution of the spatial gradients as the the big rip is approached. 
To achieve this goal Einstein equations must be written in fully inhomogeneous terms. By 
writing the line element as 2 

ds 2 = dt 2 — jij(x, t)dx l dx\ (1) 

the Hamiltonian and momentum constraints take the form 3 

K 2 — TrK 2 + r = 1QtiG[(p + p)u u° — p], (2) 
V t K - V k K* = 8nG(p + p)u° Ui , (3) 

where 

Kl = ~\ l3k Jt lk " K = Kl ^K 2 = KfK] ) (4) 
and where r = r\ is the trace of the spatial (instrinsic) curvature computed from the three 
dimensional Ricci tensor in terms of 7^. The (ij) components of Einstein equations are, 
instead, 

-L^(yyKfj - r\ = 47rG[-2(p + p)u iU i + (p - p)5{}. (5) 

where 7 = det(Ty). 

Consider then the following expansion of the spatial metric, i.e. 

7y(f,0 = a 2 {t)[aij{x) +Pij(x,t)], (6) 

when the term /3ij(x,t) contains the contribution of the gradients while cxij(x) does not con- 
tain any gradient. Recalling that the inverse metric, to this order in the gradient expansion, 

2 Note that, in this approach, jij(x,t) contains 6 independent degrees of freedom corresponding to the 
correct number of initial conditions required for a general discussion of the problem. 

3 In the following the ovcrdot will denote a partial derivation with respect to the cosmic time coordinate. 
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is given by 7* J = [a^ — P^]/a 2 (t), the extrinsic curvature (and its traces) can be readily 
computed. For instance, from Eqs. (4) and (6) 

where H = a/a. Since from the momentum constraint (3) the velocity field is always of 
higher order in the gradient expansion, i.e. 

o d k j5\ - 

U Ul = 16nG( P + P y (8) 

the contribution of the peculiar velocity field can be neglected in the remaining equations. 
This is not true necessarily to higher order in the gradient expansion. 

Thus, recalling that u u° = 1 + a lj UiUj/a 2 (t), and using Eqs. (4)-(7), Eqs. (2) and (5) 
can be written, respectively, as 

6H 2 + 2H$ + ^ = 167rGp, (9) 
a 1 

2(H + 3H 2 )5i +(3{ + 3Htt + H/35{ + \v\ = 8nG(p - p)5j, (10) 

(X 

where we defined r\ = V{ /a 2 (in this notation V\ and its traces are time-independent). 
Clearly the trace-free part of Eq. (10) reduces to 

{M - I'm) + zh(i3! - l -fj5^ = ~{n- \hv)- (ii) 

Equation (10) allows to determine the gradient contribution to the energy density and the 
following relation 

3(2H + 3H 2 )+{3 + 3Hp+^ = -24nGp, (12) 

allows to determine the gradient contribution to the pressure density. 
Even if not strictly necessary we can imagine to split p and p as 

p = p + p, p = p + p, (13) 

where, from Eqs. (10) p and p obey the usual Friedmann equations 

& = 8 -f~P, (14) 
(2H + 3H 2 ) = -8nGp, (15) 
j} + 3H(p + p) = 0. (16) 

Eq. (16) comes from the (0) component of the covariant conservation equation, i.e. 

-^4{V7[(P + P>o u ° -P\\- 4= 9 <[>/7(P + P) u o u i - K li(P + P)u k u e + p5 k e = 0. (17) 



Equation (17) also implies that p and p obey 

^ + 3H(p + p) + ^(p + p) = 0. 

According to the same logic, Eq. (9) defines p in terms of $ and V, i.e. 

2H(3 + 4 = 167rGp. 
or 

The fully inhomogeneous solution of the system can then be derived and it is: 



Ki = -H 



Si- 



te +2 



-2-3e K j 



fif(x) 



(18) 

(19) 

(20) 
(21) 



The space-dependent tensor B{ is related to the three-dimensional curvature tensor V\ by 
virtue of Eq. (10); the explicit result is 



V\ = -^[(3e + 2)(3e - 2)B\ - (3e 2 + 12e + A)B5{\. 



(22) 



The physical significance of Eq. (22) is most easily understood by inverting Eq. (22), i.e. 
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B 



#o 2 (4 - 96 2 ) 



p> - 3£2 + 12e + W 
4(9e + 4) 1 



(23) 



Equation (23) determines B\ as a function of the three-dimensional spatial curvature com- 
puted from the a^{x). The form of is in a sense arbitrary. But once «ij(f) is fixed, B\ [x) 
follow immediately from Eq. (23). 

To derive Eqs. (22) and (23) the following parametrization for the scale factor and for 
the Hubble factor 



a(t) 



to 



H(t) = H a*% 



(24) 



has been used. In Eq. (24) t br denotes the value of the cosmic time at the moment of the 
big-rip singularity. As a consequence of Eqs. (20), (21) and (22) the energy density is 



3H 2 
8^G 



l + -a 



B(x) 



(25) 



Since e > 0, the relative contribution of the gradients to the energy density (second term 
inside the squared bracket of Eq. (25)) and to the extrinsic curvature (second term inside 
the squared bracket in Eq. (21)) vanishes asymptotically, for t — > tbr- For instance, from 
Eq. (21), recalling Eq. (24) 



m = -h 



3e -\- 2 f — t 

17~ 



2(2+3e) 
3e 



(26) 



This property of big-rip singularities resembles what is known in the case of usual big-bang 
singularities where also the gradients decay though with a different power of the cosmic time 
coordinate. Big-bang singularities are argued to be homogeneous but also rather anisotropic 
[6]. As far as the big-rips are concerned, the situation may be very different. 

To illustrate the last point, consider a given time t* at which gradients can be already 
neglected. This implies, in the notation of Eq. (6), that 

qj - e 2 ^W qi - pZAyit) _ 2A z (t) (27) 

where, for simplicity, the tensor has been diagonalized 4 . 

By solving Einstein equations in the absence of spatial curvature it is easy to show that 



Mt) 



2a(U)Aj{U) 
3(e + 2) 



1 - 



-|(6+2) 



(28) 



where Ai are integration constants obeying J2iAi = 0. Initial conditions have been fixed by 
requiring that Ai(t*) = 0. Clearly, recalling Eq. (24), in the limit t — > % T the second term 
in the squared bracket of Eq. (28) goes to zero and becomes then negligible as the big-rip is 
approached. In the case of conventional big-bang singularities, the analog of the second term 
in the squared brackets evolves as a 3(w_1 )/ 2 with < w < 1. In this case the contribution of 
the anisotropy clearly grows in the limit a — > 0. 

As in the case of conventional big-bang singularities, the evolution of the anisotropy can 
also be studied in the more general case when spatial curvature is included [6]. Within 
our parametrization the contribution of the curvature always decay as a~ 2 . Therefore, for 
t — > tbr the role of the curvature may be neglected. Again, this property is not realized 
in the vicinity of big-bang singularities. Actually, in the case of big-bang singularities the 
anisotropy grows as the bang is approached. For more general classes of Bianchi models, also 
the spatial curvature grows and this occurrence may induce the celebrated chaotic features 
and the related BKL oscillations [6]. 

Up to now we have considered the case where the big-rip is caused by a perfect barotropic 
fluid with w < — 1. Needless to say that, in the class of models previously investigated, the 
dominant energy condition is violated. In connection with this problem, both in brane-world 
models [9, 10] and in four-dimensional Friedmann- Robertson- Walker models [11, 12, 13] 
the possible occurrence of a different type of quiescent (or sudden) singularities has been 
emphasized. In sudden (future) singularities the scale factor and its first derivative are 
both finite but higher derivatives of the scale factor may diverge 5 . In some class of four- 
dimensional examples this behaviour implies that while the energy density is finite at the 
rip, the pressure density diverges. As correctly pointed out in [11], sudden singularities may 
occur without a violation of the dominant energy condition. 

4 Non-diagonal (or more general) forms of ay do not change the essence of the argument (see below). 
5 Along a similar perspective but with a different chain of arguments, Ref. [14] also argues that the 
presence of a singularity in the future is not necessary even if the barotropic index is supcrnegative. 
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Consider then the example given in [11] of sudden singularities. The scale factor can be 
written as n 

a(t)=(-)% s -l) + l-(l-^) n - (29) 

For 1 < n < 2 and < q < 1 the solution is defined in the interval < t < t s . For t — > the 
model has a curvature (big-bang) singularity where the Hubble parameter and the energy 
density are both divergent. More interesting is the second singularity taking place for t — > t s . 
For t — > t B the energy density and the Hubble parameter are both finite but the pressure 
density and the second time derivative of the scale factor are divergent. 

The question we ought to address in the following concerns the nature of the gradient 
expansion in the vicinity of sudden singularities. From Eqs. (9) and (10) it can be argued 
that the contribution of the gradients to the quasi-isotropic solution may be different if the 
scale factor behaves as in Eq. (29): in the limit t — > t s both the scale factor and the Hubble 
rate are finite for the class of solutions given in Eq. (29). On the contrary, as previously 
discussed for big-rips with w < —1, in the limit t — > t\> r , the scale factor and the Hubble rate 
are divergent (see Eq. (24)). This difference is reflected in the contribution of the first-order 
gradients. 

To avoid the proliferation of parameters, the attention will now be focussed on a particular 
model belonging to the class defined by Eq. (29). Consider then the case q — 1/2 and 
n = 3/2. Indeed, the same qualitative results hold for all the models of the class described 
by Eq. (29) with the appropriate restrictions mentioned above. Let us also define, for 
notational convenience the dimensionless variable r = t/t s . Then a particular solution of 
the fully inhomogeneous system including gradients can be easily obtained and expanded for 
t ~ t s , i.e. for r ~ 1. The most notable difference is that, in this case, pressure and energy 
density are not connected by a barotropic index. The final result for the extrinsic curvature 
can be written as 

Ki = -H{t)8i-\{t)B{{x), (30) 

whith V? = t~ 2 B] and where 



H(t) = 



1 b s + 3 y/(l-T) T 



2t s 
1 



4t s a, 



1_(1_ T )2 +b s ^\ y/f 

-[2b s a s + 6a sV / T^+ (a s + 6 B )6 B (1 - r)] + 0(\1 - rf 2 ). 



6 VT^T - 3 (-5 + 4 VT^T) t + 10 b s tI + 6 t 2 v / T^7 

m = — - J , (31) 

15 f B (l - ^/T^T + b s v^+ ^/T^Tt) 



having defined, for notational convenience, b s — a s — 1. As it can be appreciated from the 
last equations, the exact expressions are rather cumbersome, therefore, in the following, the 
result for the expansion in the limit r ~ 1 will be given directly. In particular, factorizing 
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H(t) in Eq. (30) and expanding the relative contribution of the gradients we have 

I 1 [ 3 6 s a s 2 + blal V 



, (a s + 2)[l + a s 2 (4a s -13)-10a s ] 
+ 3W (1 " T) 



|3/2\ 



fi?(5) 



(32) 



Clearly, in the limit r — > 1 (i.e. t — > £ s ) the contribution of the gradients, weighted by the 
space-dependent factor B(x), is not subleading. This behaviour has to be contrasted with 
the case of big rip singularities (see, for instance, Eqs. (21) and (26)) where for t — > tbr the 
gradient contribution is subleading. 

In similar the form of p and p can also be derived and it is: 



P = P 1 



2 (a 2 - 2a s 



9 6 s 2 a s 



+ 



4 (1 + a s + a 2 



2 "2 

s "'s 

2 (-2 + 21a s + 17a s 2 + 47a s 3 -3a/ + a s 5 ) 3/2 , 

H ~ ~ Z (1 - T) + U{\1 - T\ , 



36? fl 2 



vr 



i + 



9a, 3 6 s 4 
v 7 !^ a s 2 - 1 



B(x) 



3a s 



18a 2 



l_ r ) + 0(|l_ r |3/2) 



B(x) . 



(33) 
(34) 



The exact expressions of p and p are a bit involved and then their expansion for r — > 1 will 
be given 



327rGa2 1 2 
3 



6 S 2 + 66 s v/r^ + 2Q s-5a s 2 + 13a s -l (i _ r) + _ r|3/2) 



167rGa s t s 2 v 7 !^ 



1 + 



2a? 



,(35) 
(36) 



From Eq. (33) it can be deduced that as t — > t s the gradients are not subleading. Moreover, 
from Eq. (35) one can also argue that for t — > t s , p is finite. On the contrary, in the same 
limit, i.e. r — > 1 the pressure density diverges (see Eq. (36)). The amusing thing is that, in 
this case, the relative contribution of the gradients of Eq. (34) is subleading as t — > t s . 

In conclusion, let us summarize the main findings of the present investigation. In the 
first part it has been shown that if the dominant source of the background geometry is given 
by a perfect fluid with supernegative barotropic index (i.e. w < —1) then the contribution 
both of the spatial gradients and of the anisotropy tends to decay as the big-rip singularity is 
approached. For the more conventional big-bang singularities the situation is a bit different: 
while gradients also decay in the vicinity of the big-bang, the anisotropy and the curvature 
may well grow and lead to some type of chaotic approach to the singularity. 

We then moved to the analysis of sudden singularities. In this case the dominant energy 
condition is not violated. While the scale factor, the Hubble parameter and the energy 
density are all finite as the singularity is approached, the pressure density diverges. In this 
situation we included the contribution of the gradients and showed in an explicit example 
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(representative of a more general class of backgrounds) that the relative contribution of the 
gradients does not decay in the vicinity of the sudden singularity. 

Various interesting generalizations are left for future works. The analysis of single fluid 
big rip singularities can be generalized to a multi-fluid situation in analogy to what recently 
discussed in the homogeneous and isotropic case [15]. In this situation it would be also 
interesting to discuss the fate of the anisotropy in more general Bianchi models. 
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